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Abstract 

We suggest theoretical explanation of the high upper critical magnetic field, perpendicular to conducting chains, 
H c2 , experimentally observed in the superconductor (TMTSF) 2 C104, in terms of singlet superconducting pairing. 
In particular, we compare the results of d-wave-like nodal, d-wave-like node-less, and s-wave scenarios of super- 
conductivity. We show that, in d-wave-like nodal scenario, superconductivity can naturally exceed both the orbital 
upper critical magnetic field and Clogston-Shandrasekhar paramagnetic limit as well as reach experimental value, 
H C 2 — 6 T, in contrast to d-wave-like node-less and s-wave scenarios. In our opinion, the obtained results are strongly 
in favor of d-wave-like nodal superconductivity in (TMTSF)2C1C>4, whereas, in a sister compound, (TMTSF)2PF6, 
we expect either the existence of triplet order parameter or the coexistence of triplet and singlet order parameters. 
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1. Introduction 

High magnetic field properties of (TMTSF) 2 X 
(X=C104, PF6, etc.) organic materials have been 
intensively studied [1] since the discovery of super- 
conductivity in (TMTSF) 2 PF 6 conductor [2]. From 
the beginning, it was clear that superconductivity in 
the above mentioned materials was unconventional. 
Indeed, early experiments [3,4] showed that the Hebel- 
Slichter peak was absent in the NNR measurements 
[3] and that superconductivity is destroyed by non- 
magnetic impurities [4] . These facts were strong argu- 
ments that superconducting order parameter changed 
its sign on the quasi-one-dimensional (Q1D) Fermi 
surfaces (FS) of (TMTSF) 2 X material. The main re- 
sults of both experiments [3,4] were recently confirmed 
in a number of publications (see, for example, [5,6]). 
It is important that the above mentioned experiments 
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did not contain information about spin part of super- 
conducting order parameter and could not distinguish 
between singlet and triplet superconducting pairings. 

The first measurements of the Knight shift in 
(TMTSF) 2 PF 6 conductor [6,7] showed that it was not 
changed in superconducting phase, which was inter- 
preted in favor of triplet superconductivity [6,7]. On 
the other hand, the more recent Knight shift mea- 
surements [8] in superconductor (TMTSF) 2 C1C>4 have 
shown a clear change of the Knight shift in super- 
conducting phase at relatively low magnetic fields, 
H ~ 1 T, and have been interpreted as evidence of sin- 
glet superconductivity [8]. Another argument in favor 
of singlet order parameter is the fact that the upper 
critical magnetic field, parallel to conducting axis, H^ 2 
[9], is paramagnetically limited [10]. Moreover, very 
recently the Larkin-Ovchinnikov-Fulde-Ferrell phase 
[11,12], which appears for singlet superconducting 
pairing, has been experimentally discovered [13,14] in 
(TMTSF) 2 C10 4 and theoretically interpreted [15]. 



Preprint submitted to Physica B 



3 September 2012 



2. Goal 

In such situation, where support of singlet super- 
conducting pairing in (TMTSF^CICU material is in- 
creasing, it is important to reinvestigate theoretically 
high experimental upper critical magnetic fields, H^2 
[13,14,16], observed for a magnetic field, perpendicular 
to conducting chains. For many years, large values of 
have been considered as a consequence of triplet 
superconducting pairing. Our goal is to show that we 
can naturally explain large values of H^ 2 within sin- 
glet d- wave- like nodal scenario of superconductivity in 
(TMTSF) 2 C10 4 . We also show that d-wave nod-less 
and s-wave scenarios are much less consistent with ex- 
perimental value of Hc2 - 67 [13,14]. This value ex- 
ceeds both the quasi-classical upper critical field [17] 
and Clogston-Shandrasekhar paramagnetic limit [18] 
due to the coexistence of two unusual superconduct- 
ing phases: Reentrant superconductivity [19-24] and 
Larkin-Ovchinnikov-Fulde-Ferrell phase [11,12]. 



3. Results 

Let us consider Q1D spectrum of (TMTSF) 2 C10 4 
conductor in tight binding model [1], 

e(p) = — 2t a cos(p x a/2) — 2t b cos(p y b') — 2t c cos(p z c*),(l) 

in a magnetic field, perpendicular to its conducting 
chains, 



H = (0,/J,0), A = (0,0, —Hx), 



(2) 



where t a ^> tb^> t c are electron hoping integrals along 
a , b , and c* axes, respectively. Electron spectrum (1) 
can be linearized near two slightly corrugated sheets of 
Q1D FS as 

^(p) = ±Vx(p y )\px T Pf(Pv)] - 2t c cos(p z c*), (3) 

where +(-) stands for right (left) sheet of Q1D FS. 

We represent electron wave functions in a real space 
in the following way: 

yf(x,y,z,a) = exp[ip F (p y )x] exp(ip y y) exp(ip z z) 

Xipt(x,p y ,p z ,a). (4) 

Let us use the Peierls substitution method, 

Px TPf(p v ) -»■ -id/dx, p z -»• p z - eA z /c. (5) 

As a result, the Schrodinger-like equation for wave 
functions ^(x, p y , p z , a) can be written as 



^iv x (p y )-^ 2t c cosl p z c* + —x ]—Li B aH 

ax \ vf 



Xtp e (x,p y ,p z ,a) = 8e ip e (p x ,p y ,p z ,(ff) 



where /is is the Bohr magneton, a = ±1 stands for 
spin up and down, respectively; lu c = eHv F c* /c, 8e — 
a - e F . 

It is important that Eq.(6) can be analytically 
solved: 



i>e (x,p y ,p z ,a) = 



exp[±i5ex/v x {p y )\ 



y/2-KV x {p y 



exp 



±i 



\iboHx 



V X (Py 



x exp 



±i 



2t c 



V x (Py) 



|cos(; 



p z c* H -u du 

v F 



(7) 



The corresponding Green functions can be obtained 
from the following equation (see Ref.[25]): 



G(x,x 1 ,py,p z ,a) = 



ip*(x,Py,Pz,cr)ipe{xi,Py,p z ,a) 



.(8) 



Below, we introduce superconducting order param- 
eter in the following way: 

A(p y ,x) = f(p y b')A(x), J f(p y b')d(pyb')/2iT = 1,(9) 



where function f(p v b') takes into account three pos- 
sible order parameters: d- wave-like nodal, f(p v b') = 
y/2 cos(p y b'), s-wave, f(p y b') = 1, and d-wave- like 
node- less, f(p y b') = 28(p y b' + n/2) - 26{p y b' - n/2) - 1 
[26] . Linearized gap equation for all three possible sin- 
glet scenarios of superconductivity in (TMTSF)2C104 
conductor can be derived, using the Gor'kov equations 
[25] for non-uniform superconductivity [28-30]. As a 
result of rather lengthy but straightforward calcula- 
tions, we obtain: 



J V'APy) J 



2nTdxi 



, , 8t c V F 
X JoS ; 7 sm 

tOcV x (Py) 



X — Xl | > jy^— 

UJ c {x — Xl) 



^(P,)sinh[2^i 



2v F 



LU c (x + Xl) 



2v F 



X cos 



2/3hbH(x - xi) 



V X (Py) 



f{ P yb')A{Xl) 



(10) 



where g is effective electron coupling constant, Q. is cut- 
off energy, parameter /3 takes into account possible de- 
viations of superconductivity in (TMTSF)2C104 from 
weak coupling scenario. 

Note that Eq.(10), derived in this article, is very gen- 
eral. It takes into account both the orbital and para- 
magnetic effects against superconductivity. In partic- 
ular, it takes into account quantum nature of electron 
motion along open orbits in the extended Brillouin 
zone in Q1D conductor (3) in a magnetic field (2). It 
is possible to make sure that the main quantum pa- 
rameter in Eq.(10) is 2t c VF /w c v x {py) — 2t c /u) c - Let us 
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estimate the value of this quantum parameter, using 
quasi-classical language. In accordance with Ref.[19], 
quasi-classical electron trajectory in a magnetic field 
(2) can be written as 



z(t,H) = c l±(H)cos(u c t), 



(11) 



where l±(H) = 2t c /u) c corresponds to a "size" of elec- 
tron trajectory in terms of interlayer distance, c*; t is 
time. 

It is possible to show that 



2V2 4>o t c ^ 2 x 10 3 t c h 
~iT ac*HU ~ H(T) 



(12) 



where H(T) is a magnetic field, measured in Teslas, 
4>o is the flux quantum. Note that value of t a /tb — 10 
is very well known in (TMTSF^CICU from theoretical 
fitting [31] of the so-called Lee-Naughton-Lebed oscil- 
lations [31,32]. As to ratio t c /tb, it can be evaluated 
from the measured Ginzburg-Landau (GL) slopes of 
the upper critical magnetic fields in (TMTSF)2C104 
conductor [13,14]: 



t c /t b = (b*/V2c*-)(H c c2 /H b C 2) G L, 
tc/h = (b*/c*)(H c c2 /H%) G L, 



(13) 
(14) 



where Eq.(13) is valid for d- wave- like nodal pairing, 
whereas Eq.(14) is valid for both d- wave-like node- less 
and s-wave pairings. 
As a result, we obtain 

l±(H = 6 T) ~ 0.48 (15) 

for d- wave- like nodal scenario of superconductivity and 



l±(H = 6 T) ~ 0.68 



(16) 



for d-wave-like node-less and s-wave ones. Low values 
of the parameter l±(H = 6 T) show that both cases 
correspond to 3D — > 2D dimensional crossover of elec- 
tron motion [19], where electrons are almost localized 
on conducting layers. In this situation, superconduc- 
tivity becomes almost two-dimensional and we can ap- 
proximate the Bessel function in Eq.(10) as Jo (2) ^ 
2 2 /4. 

Below we consider the gap equation (10) at zero tem- 
perature, T = 0. In this case and under condition of 
3D — > 2D dimensional crossover, it is possible to repre- 
sent superconducting order parameter in the following 

way: 



A(x) = exp(ikx)[l +ai cos(2uJcx/vf) 
+a 2 sin(2w c a;/i>F)], 



(17) 



where Qi|, Q2I <SC 1. For such order parameter Eq.( 10) 
can be rewritten as 



1 
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Vf 



2tt/6' 



/ d(p y b') f dz 2 , / 

J J T /feb)cos 



2/3ii b Hz 
y v F 



Vxijpy 

v F 



CIS) 



where g is renormalized electron coupling constant, 
xi—x — zv x (p v )/vf- Here, we consider Eq.(18), taking 
into account that electron velocity component along 
conducting x axis is 



Vx(p y ) =v F [l + acos(p y b')], 



(19) 



where a = V2t b /t a ~ 0.14 [23]. Under condition a <ti 
1, Eq.(18) can be simplified: 



00 

1 _ r dz 

9 J z 



L-2Ji(ff)Bin 2 (fg| , m 



where m = 1 for d-wave-like nodal superconducting 
order parameter, whereas m — for d-wave-like node- 
less and s-wave ones. It is important that, in the ab- 
sence of the paramagnetic effects, Eq.(20) describes the 
Reentrant superconductivity [19] with transition tem- 
perature being increasing function of a magnetic field. 
Therefore, we call superconducting phase, described by 
Eq.(20), the hidden Reentrant superconductivity. 

We can simplify Eq.(20) by using the following rela- 
tionship: 



oc 

9 J 



2nT c dz 



vf sinh(- 



(21) 



where T c is superconducting transition temperature in 
the absence of a magnetic field. As a result, we obtain 



In 



H b c2 \ f dz I 

IF =J T cos { 
' x 



2/3/ibHz 



vf 



cos(kz){Jo(akz) — mJ 2 (akz)\ 



2/i(7J)sin 2 



2v~f~ J 



(22) 



where [ibH* = nT c /2 r y, 7 is the Euler constant. We 
numerically find maxima of H b 2 (k) as a function of 
wave vector, k, of the order parameter (17) for m=l 
and m=0. We come to the conclusion that experimen- 
tal value, H h c2 ~ 6 T, can be obtained for d-wave-like 
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nodal order parameter at /3 ~ 0.85, whereas for d- 
wave-like node-less and s-wave order parameters it cor- 
responds /3 ~ 0.5. Note that the so-called (/-factor in 
(TMTSF)2C104 conductor is very close to its standard 
value, g= 2 [1], which corresponds to j3 = 1 in Eq.(22). 
Therefore, we consider d- wave-like nodal superconduc- 
tivity to be much more consistent with the experimen- 
tal data, where the calculated value /3 ~ 0.85 is close 
to its expected value, j3 — 1. It is not exactly equal to 
1, perhaps, due to slightly deviations of superconduc- 
tivity in (TMTSF)2C104 from weak coupling scenario. 
On the other hand, we consider d-wave-like node-less 
and s-wave parings as very unlikely due to very low 
calculated value of /3 ~ 0.5. 



4. Conclusion 

We have suggested explanation of high values of 
the upper critical magnetic fields, experimentally ob- 
served in (TMTSF) 2 C10 4 conductor [13,14,16], using 
d-wave-like nodal scenario of superconductivity. On 
the other hand, we anticipate that, for explanation of 
very high values of the upper critical magnetic fields 
[33,34] in a sister compound, (TMTSF) 2 PF 6 , in mixed 
superconducting-antiferromagnetic phase [34,35], it is 
necessary to consider either triplet superconducting 
pairing [10,23,24] or the coexistence of triplet and sin- 
glet superconducting order parameters [36] . 

We are thankful to N.N. Bagmet, S.E. Brown, and Y. 
Maeno for useful discussions. This work was supported 
by the NSF under Grants No DMR-0705986 and DMR- 
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